Lukewarm black holes in quadratic gravity. 
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Perturbative solutions to the fourth-order gravity describing spherically-symmetric, static and 
electrically charged black hole in an asymptotically de Sitter universe is constructed and discussed. 
Special emphasis is put on the lukewarm configurations, in which the temperature of the event 
horizon equals the temperature of the cosmological horizon. 
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I. INTRODUCTION 

o ■ 

fNl ' A complete description of the gravitational phenomena should be given by the quantum gravity, which, in turn, 
may be a part of the even more fundamental theory. Unfortunately, at the present stage we have no clear idea how 
^ ■ such a theory should be constructed. It is expected, however, that the action functional describing its low-energy 
approximation should consist of the higher order terms constructed from the curvature and its covariant derivatives 
to some required order. Such generalizations of the Einstein-Hilbert action functional necessarily introduce a number 
of coupling constants, which, at the present state of affairs, should be determined empirically. 

For traces of possible breakdown of the classical General Relativity one, quite naturally, turns to cosmology. Indeed, 
it is possible that such concepts as dark matter or/and dark energy proposed in order to explain observational data 
should be abandoned in favour of modifications of the classical gravitational Lagrangian (see e.g. Ref. 1 and references 
^ ' cited therein). 

(3J[)! Among various modifications of the general relativity proposed so far, a prominent role is played by the quadratic 
>— gravity. Various aspects of such theories have been discussed extensively in the literature. (See for example Refs. 
2-12). As is well-known the motivations for introducing linear combination of R 2 , R a bR ab and R a bcdR ahcd into the 
gravitational action functional are numerous. For example, when invented, the equations of quadratic gravity have 
■ been treated as an exact formulation of the theory of gravitation. For historic informations and important references 
the reader is referred to Ref. 13. It may be considered, quite naturally, as truncation of series expansion of the action 
• of the more general theory. Moreover, from the point of view of the semi-classical gravity it might be treated as some 
sort of a poor man's stress- energy tensor, allowing in a relatively simple way to mimic, especially when the application 
\Q • of the full stress- energy tensor would produce extremely complicated result, the fairly more complex source term 
of the field equations. This is especially true when the right-hand-side of the semiclassical Einstein field equations 
is taken to be the renormalized stress-energy tensor of the quantized fields in a large mass limit constructed within 
the Schwinger-DeWitt framework. Such tensors comprise the linear combination of the purely geometric terms and 
the spin of the field enters through the numeric coefficients fl4l - fl8| . As compared with the General Relativity the 
Lagrangian of the quadratic gravity in four dimensions requires two additional terms aR a bR ab + PR 2 , where a and 
j3 are the coupling constants. The possible third term constructed form the Kretschmann scalar can be relegated by 
means of the Gauss-Bonnet term. 

In general, the equations derived form the higher-order action functional are very hard to solve. Fortunately, since 
the coupling constants are expected to be small one can easily employ a perturbative approach to the problem treating 
the classical solution of the Einstein field equations as the zeroth-order of the approximation. Successive perturbations 
are therefore solutions of the differential equations of ascending complexity. It should be noted that although the 
method is clear the calculations beyond the first- order may by intractable. Such prerturbative approach is in concord 
with the philosophy of the effective theories. 

Exact and approximate solutions to the equations of the quadratic gravity have been studied in a number of papers. 
Of particular interest are the spherically-symmetric configurations which, potentially, may describe black holes. The 
black hole solutions have been studied for various sources in Refs. 19-24. Here we shall analyze a particular class of 
solutions of the equations of the quadratic gravity with the cosmological term describing lukewarm black holes [25142 71 ]. 

For the Reissner-Nordstrom-de Sitter class of solutions it is possible to find configurations in which the event and 
cosmological horizons have the same temperature. These solutions are known as lukewarm black holes and the interest 
in them stems from the fact that if Th = T c it is possible to construct a regular thermal state for a two dimensional 
models [28| . Moreover, analyses of the field fluctuation indicate that this is the case in four dimensions [291 [30j] . (It 
should be noted that such configurations are prohibited in the geometry of the Schwarzschild-deSitter black hole.) The 
natural question that arises in connection with the foregoing discussion is whether or not it is possible to construct 
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the lukewarm black hole in the quadratic gravity. And although the full, detailed answer is beyond our capabilities, 
it is possible to provide an affirmative answer to the restricted problem. Indeed, since the complexity of the coupled 
equations of the quadratic gravity and electrodynamics, even in the simplest case of spherically-symmetric and static 
geometries, hinders construction of the exact solution, one has to refer to the analytical approximations or numerical 
methods. Here we shall employ the perturbative approach. The obtained results can also be viewed as a first step 
towards incorporation of the quantum effects into the picture. It is because the renormalized stress-energy tensor of 
the quantized massive field may be approximated by the object constructed from the curvature tensor, its covariant 
derivatives and contractions. 



II. EQUATIONS OF THE QUADRATIC GRAVITY 

The coupled system of the classical electrodynamics and the quadratic gravity is described by the action 

S = ^S g + S m (1) 



where 



and 



S g = J (R + 2A + aR 2 + (3R ab R ab ) \[^~g d 4 x, (2) 



S m = -T^- / F ab F ab V=gd A X, (3) 
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where all symbols have their usual meaning. The third possible term constructed form the Kretschmann scalar, 
RabcdR abcd , may by removed from the Lagrangian with the help of the Gauss-Bonnet invariant 

RabcdR abcd -4R ab R ab + R 2 . (4) 

Of numerical parameters a and (3 we assume, as usual, that they are small and of comparable order, otherwise they 
would lead to the observational consequences within our solar system. Their ultimate values should be determined 
form observations of light deflection, binary pulsars and cosmological data [3~iT-l33j| . Moreover, following Ref. 19, we 
shall restrict ourselves to spacetimes of small curvatures, for which the conditions 

\aR\ < 1, \pR ah \ < 1 (5) 

hold. Although the additional constraints could be obtained from non-tachyon conditions [l2l. l2lj|. which, in turn, are 
simple consequences of demanding the linearized equations to posses a real mass, we shall treat the parameters a and 
f} as small but arbitrary. 

Differentiating functionally the action S with respect to the metric tensor one has 

L ab = G ab + Agab _ a jab _ pjab = g^yab^ (g) 

where 

jab = 2R ; ab _ 2RR ab + }_ g ab ( R 2 _ 4Q ^ ^ 

and 

jab = R - ab _ aR ab _ 2R cd R chda + ^g ab (R cd R cd - UR) . (8) 

Let us consider the spherically symmetric and static configuration described by the line element of the form 

rlr 2 

ds 2 = _ e mr) f{r)dt 2 + |!_ + r W, (9) 

fir) 



where 
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The spherical symmetry places restrictions on the components of F a b tensor and its only nonvanishing components 
compatible with the assumed symmetry are i*oi an d F23. Simple calculations yield the stress-energy tensor in the 
form 

mt _ rpr _ Q 2 + P 2 _ Z 2 

and 



n 2 4- p 2 

T^T^^ S ^, (12, 



where the integration constants Q and P are interpreted as the electric and magnetic charge, respectively. 

III. PERTURB ATI VE SOLUTION 
A. General Case 

To simplify calculations and to keep control of the order of terms in complicated series expansions we shall introduce 
a dimensionless parameter e substituting a — > ea and f3 — > ef3. We shall put e = 1 at the final stage of calculations. 
Of functions M (r) and ip (r) we assume that they can be expanded as 

M (r) = M (r) + eMi (r) + O (e 2 ) (13) 

and 

i> (r) = &h. (r) + O (e 2 ) . (14) 

Consider the left hand side of Eq. ^ calculated for the line element © first. Making use of the above expansions 
and collecting the terms with the like power one obtains 
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where 



L* = A-^(Mj+eMj-eSf), (15) 



S* = p [ EM _ 8M o M o + 2M o 2 _ + 5M Aff _ M' Q M^ 

(3 ) _ MoMf _ M , ( 3 ) + r M ( 4) _ 2 Mq M (4)\ 



M" 2 MnM n (3) 
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(2AM0K 8M 4M; 2 8Af,? 18 M MX „ a 



- 4 Mf + 6M ° r M ° 3) + 2 Mq M< 8 > - 2 r M< 4 > + 4 M M< 4 ^ (16) 

(i) 

and primes as well as Mq for i > 3 denote derivatives with respect to the radial coordinate. It should be noted that 
Eq. (T151) does not depend on the function ip(r). The zeroth-order and the first-order equations integrated with the 
conditions Mo(r + ) = r + /2 and Mi(r + ) = 0, respectively, give 

m = 1 - - 

r 



(17) 
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It should be noted that a correct choice of r + in the zeroth-order solution requires some prescience. Indeed, for a given 
Z 2 and A it should allow for a positive root, say c, interpreted as the cosmological horizon which satisfies r + < c. 
Remaining roots, say a and b, should satisfy a < b < r + < c. 

On the other hand, the difference between the radial and time component of the tensor L b a can be easily integrated 
to yield 

Mr) = (2a + /3)a4 3) ~^(3a + $)M' + C x , (18) 
where C\ is the integration constant. To determine C\ we shall adopt the natural condition [34j 

9tt (roc) 9rr (roo) = -1, (19) 

or, equivalently, 

^(roo) = 0, (20) 

where r m is either infinity or the cosmological horizon or the radius of a cavity in which the system is placed (the 
latter case is not considered in this paper). Simple integration gives 

and the equations (|17l) and (1211) provide complete first-order solution to the problem. In what follows, for simplicity, 
we shall choose r m = oo. The other choices do not influence location of r c nor the relations describing equality of 
temperatures of the cosmological and the event horizons. 



The zeroth-order solution can be expressed in a more familiar form using, for example, the Abbott-Deser mass 
Indeed, if the coefficient that stands in front of r _1 is small it is possible to relate it with the Abbott-Deser mass 
On the other hand, one can introduce the parameter 



to get 



Simple manipulations in Eq. (|22[) give 



M = ^ + ^-M (22) 
2 2r+ 6 v ' 



/w = i ~—+%- A?- m 

r r 6 



2M Q 2 AH 



r+ r\ 



(24) 



and the parameter M. can be referred to as the horizon defined mass. 

Returning to the first order solution one concludes that the cosmological horizon is given by 



£ r«, (25) 



where 



(1) 4AQ 2 (f c -r + ) 3(f c 5 -5r4f e + 4r^)Q 4 

' V - •■ a "i i ,, • i. .-, P 



(-4 f c r\ + A f A c r\ - A r\ + 3 r\ + ff) Q 2 



(26) 



and k c = /'(f c )/2, provided f c is the greatest positive root of f(r) — 0. The function f(r) is the unperturbed part of 
f(r). Equations (fTTT) . (|2T|) and ([26)1 constitute the full (first-order) solution solution to the problem. 
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B. T H = T c case 



Now, we shall restrict ourselves to the particular class of solutions describing black holes for which temperature 
of the event horizon equals the temperature of the cosmological horizon. Consider the zeroth-order solution (e = 0) 
first. Solving the system 



/(f c ) = and f'(r p ) + f'(r c ) = 



with respect to A and Z 2 one has 



A 



r~) 2 



and Z 2 = 



r+r c 



J+ + r c 

In the (r + , r c )— parametrization the line element assumes the form (j9|) with ip(r) = and 

\ 2 2 

r + r c x 



m = (i 

whereas in the (I, z) — parametrization one has 



(r+ + f c )r J (r+ + f c ) 2 



1 



r 



where z = \/Q 2 + P 2 and Z = y/3/A. For 4z < Z one has two pair of roots: 

and f r — 
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which are interpreted as the event and cosmological horizons, and 
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(32) 



The positive root r_ is interpreted as the inner horizon whereas r is a negative root which has no physical meaning. 

By the construction the temperature of the event horizon equals that of the cosmological horizon 



T H = T C = 



2ir(r c 



V2 ' 



(33) 



Such configurations are usually referred to as the lukewarm black holes [2514271 l29l |30| • From the point of view of 
the quantum field theory in curved background the lukewarm black holes are special. It has been shown that for the 
two-dimensional models it is possible to construct a regular thermal state [28| . Moreover, recent calculations of the 
vacuum polarization indicate that it is regular on both the event and cosmological horizons of the D— 4 lukewarm 
RN-dS black holes gHH. 

The lukewarm configuration in the quadratic gravity should simultaneously satisfy 



f(r c ) = and T H = T c . 



(34) 



The Euclidean section is regular at r+ and r c if r = it is periodic with a period 2tt /k, where k is the surface gravity. 
The equality of the temperatures is equivalent to 



1 



\r=r+ 



1 



d 



0. 



(35) 



and, consequently, the lukewarm configuration is described by the above equation and the first equation of (|34|) . Since 
r+ always denotes the exact location of the event horizon, mathematically, we have two relations for three parameters 
and to describe the lukewarm configuration completely Eqs. (|34l) should be supplemented by additional condition. 
Here, following Ref. 32, we shall adopt the point of view that the cosmological constant is not a parameter of the 
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space of solutions, but, rather, the parameter in the space of theories. To construct the perturbed lukewarm black 
hole let us substitute (123 and 



and 



^ = UrU +£A (36) 

into the line element. Solving the system one obtains 

2A? 2 c r\a 2r + (5 r\ + 20 rfrj. + 28 r%r+ + 6 f c Hj. + r\) /3 
[r+ + r c ) 5 (r+ + r c J 

_(i) _ g (^-r + ) 2 (fg+47V + +r^) 
or c [r+ +r c ) r + 

The higher-order curvature corrections to the geometry of the classical lukewarm black holes can readily be obtained 
by substituting (|37[) and ([38} into the line element ((9]), and, after expansion, retaining the terms linear in s. 

It is of some interest to examine a few special cases of the results (|36M38p . In the limit r c — > oo the correction 
A tends to zero and one obtains the first order corrections to the extreme Reissner-Nordstrom solution with the 
near-horizon geometry of the type AdS2 X S 2 . If j3 — and a^O the correction to the cosmological horizon vanishes 
and 

9 ( r+r c V 24f>ia 

Z = + , -n4 - ( 39 ) 

\r++r c J (r++r c ) 

On the other hand, if a = and (3^0 both A and r c are always of the same sign as /?. Finally, if 3a + (3 — both A 
and r c are always of the same sign as —a. Additional constraints can be obtained from the non-tachyon conditions. 



IV. FINAL REMARKS 



In this paper we have constructed the first order-correction to the geometry of the Reissner-Nordstrom-deSitter black 
holes within the framework of the quadratic gravity. Special emphasis has been put on the lukewarm configurations 
which are characterized by the equality of the temperature of the event and cosmological horizon. Though we have 
studied the consequences of the inclusion of the simplest higher-order correction to the gravitational action our 
results can shed some light on the important issue of the semi-classical lukewarm black holes. Since the approximate 
renormalized effective action of the quantized massive fields in the large mass limit is constructed solely form the 
curvature invariants (the type of the field is encoded in the numerical coefficients) the analysis should follow the same 
lines as the analysis presented here. This strongly support the hypothesis that solutions describing the quantum- 
corrected lukewarm black holes do exist. We intend to return to this group of problems elsewhere. 
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